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ABSTRACT 

We present the results of a three years monitoring campaigns of the z = 0.024 type-1 active galactic nucleus (AGN) PGC50427. Using 
robotic telescopes of the Universitatssternwarte Bochum near Cerro Armazones in Chile, we monitored PGC50427 in the optical and 
near-infrared (NIR). Through the use of Photometric Reverberation Mapping with broad and narrow band filters, we determine the 
size of the broad-line emitting region by measuring the time delay between the variability of the continuum and the Her emission 
line. The Her emission line responds to blue continuum variations with an average rest frame lag of 19.0 + 1.23 days. Using single 
epoch spectroscopy obtained with the Southern African Large Telescope (SALT) we determined a broad-line Her velocity width of 
1020 km s"^ and in combination with the rest frame lag and adoption a geometric scaling factor / = 5.5, we calculate a black hole 
mass of Mbh ~ 17 X 10^M©. Using the flux variation gradient method, we separate the host galaxy contribution from that of the 
AGN to calculate the rest frame 51OOA luminosity at the time of our monitoring campaign. We measured small luminosity variations 
in the AGN (~ 10%) accross the 3 years of the monitoring campaign. The rest frame lag and the host-subtracted luminosity permit 
us to derive the position of PGC50427 in the BLR size - AGN luminosity diagram, which is remarkably close to the theoretically 
expected relation of R oc LP-^. The simultaneous optical and NIR (/ and K^) observations allow us to determine the size of the dust 
torus through the use of dust reverberation mapping method. We find that the hot dust emission (~ 1800X) lags the optical variations 
with an average rest frame lag of 46.2 + 2.60 days. The dust reverberation radius and the nuclear NIR luminosity permit us to derive 
the position of PGC50427 on the known r - MV diagram. The simultaneus observations for the broad-line region and dust thermal 
emission demonstrate that the innermost dust torus is located outside the BLR in PGC50427, supporting the unified scheme for AGNs. 

Key words, galaxies: active -galaxies: Seyfert -quasars: emission lines -galaxies: distances and redshifts -galaxies: individual: 
PGC50427 


1. Introduction 


The existence of a dusty structure with a torus-like geometry sur¬ 
rounding the broad-line region (BLR) and the central accreting 
super massive black hole (SMBH) play a fundamental role in the 
framework of an unified model for active galactic nuclei (AGN; 
AntonuccI||1993| ). Its presence would explain the observed dif¬ 
ferences in the spectra of type 1 and 2 Seyfert galaxies together 
with the strong emission observed at near/mid/far infrared wave- 
lengths ([Alonso-Herrero et al.|2001t [Haas et al.|2003t[Haas et aL] 


20Q4[[sTebenmorgen et al.|2065|). 

Evidence for the surrounding hot dust torus rests on indi 
rectly optical spectropolarimetric observations (eg., |Antonucci 


larvainis||1987t [Kobayashi et~ar|1993|). A di- 


rect detection of the putative dust torus is more difficult since the 
internal structure of AGNs is spatially unresolved using single 
telescopes. IR long- baseline interferometric observations have 
been able to determine the size of the dust torus for a few nearby 


* Based on observations made with the Southern African Large Tele¬ 
scope (SALT) 


AGNs (eg.,|Wittkowski et al.|2004l [Tristram et al.|2007l|Pott et 
|al.|2010D . 

The only method available to study the origin and morphol¬ 
ogy of the BLR and the dust torus independent of the spatial res¬ 
olution of the instrument is Reverberation Mapping (RM, [Lyutyi 


& Cherepashchuk|1972[|Cherepashchuk & Lyutyi|1973 [Bland- 


ford & McKee|1982t|Gaskell & Sparke|1986t[Peterson T993t|Pe- 


terson et al.|2()04[ ). In this method one measures the light travel 

time between the accretion disk (AD) and the BLR and/or hot 
dust. The hot AD produces a variable continuum emission, and 
this variability is observed with a time delay (tblr = RblrIc) in 
the studied broad emission lines of the BLR. Similarly, if the dust 
torus is located at some radial distance Rdust around the hot AD it 
will reprocess the UV/optical radiation to thermal near-infrared 
(NIR) radiation with a characteristic time delay Tdust = ^dustlc. 
Reverberation mapping has revealed the size of the BLR, the 
black hole mass and Eddington ratios in about 50 AGNs (see [Du[ 
[et al.|2014| and references therein). 

The torus reacts to the AD variability in a wide range of 
time, from days to decades (eg., [Glass[[2004t [Suganuma et al 


[2006[[Pozo Nunez et al.[[2014|[Koshida et al.|2014[ ). Erom dust 

Article number, page 1 ofpd] 










































A&A proofs: manuscript no. dust_BLR 


VYSOS-6 

VYSOS-6 

VYSOS-6 


BEST-II 

BEST-II 

IRIS 

IRIS 


BEST-II 

BEST-II 

BMT 


Fig. 1. The timeline of our campaign on PGC50427. The days in which the telescopes made observations are indicated by crosses. 
Table 1. Characteristics of PGC50427 
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Values from NED database, Veron-Cetty & Veron ( 201Q| , Schlafly & Finkbein^|201l| . 


reverberation measurements, which are principally based on the 
cross-correlation analysis between the optical (V, 0.55 pm) and 
NIR (K, 2.2 pm) light curves, a correlation between the inner¬ 
most radius of the torus and the squ are root of the optical lumi- 
nosity has been shown {Rdust ^ Glass|20Q4 Minezaki et al 


|20Q4| Suganuma et al.|[2006} Koshida et al.||2014| ). The dust is 

heated by the accretion disk (AD) until up to its maximum subli¬ 
mation temperature (~ 1500K). The correlation and sublimation 
temperature are consistent with the theoretical prediction of the 
dust sublimation radius {Rsub^ |Barvainis|199^ , but is systemat¬ 
ically smaller by a factor of three than the sublimation radius 
Rsub predicted for graphite dust grains with a size of 0.05 pm 
in radius ( [Qknyanskij et al.|1999||Kishimoto et al.|2007| ). Some 
modihed dust geometries involve BLR associated dust due to 
winds/outflows from the accretion disk ([Konigl & Kartje|1994 


Elitzur & Shlosman|2006[[Czerny & Hryniewicz|2011 ). 


Alternatively, theoretical simulations shows that an 
anisotropicahy illuminated dust torus, caused by an optically 
thick AD, places the inner concave region of the torus closer 
to the outer edge of the AD. Thereby increasing the response 
time of the torus, and thus explaining the systematic difference 
of the time delay with respect of the torus radius measured 
from time delay and from the sublimation temperature under an 
isotropically a ssumption ( [Kawaguchi & Mori||2010||Kawaguchi| 


& Mori||2Qll|). In the c ase of the Seyfert 1 galaxy WPVS48, 
Pozo Nunez et al.| ( [2014| ) argue that the sharp NIR echo observed 


is due to a geometrically and optically thick torus seen nearly 
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face-on. In this scenario the observer only sees the facing 
rim of the torus wall, which lies closer to the observer than 
the torus equatorial plane and therefore leads to an observed 
foreshortened lag effect. 

While great theoretical progress has been made, only a hand¬ 
ful of observational measurements of the dust reverberation ra¬ 
dius have been obtained during the last years. Moreover, as noted 
in [Suganuma et al.| ( |2006| ), the importance of simultaneous BLR 
and dust torus size measurements provide an important step for¬ 
ward to test and constrain the actual paradigm of unihcation for 
AGNs. 

More recently. Photometric Reverberation Mapping (PRM) 
has been revisited and used as an efficient alternative to de¬ 
termine the BLR s ize, black hole masses and host-subtracted 
AGN luminosities ([Haas et al.[[2Qll[ [Pozo Nunez et al.[[2Q12[ 
[Pozo Nunez et al.|2013 ). Through the combination of broad and 
narrow-band data, this method is used to measure the time de¬ 
lay between the triggering continuum variations and the BLR 
emission line response, which has previously been isolated by 
the subtraction of the underlying continuum determined from the 
broad-band hlter data and/or through a single spectrum contem¬ 
poraneous with the campaign. 

Feature-rich PRM light curves allow us to infer the basic ge¬ 
ometry of the BLR, whether it is spherical or disk-like, and can 
thus constrain the unknown geometrical factor needed in con¬ 
verting the time lag and velocity width into a black hole mass 
( [Pozo Nunez et aL][2014[ ). Furthermore, the use of broad-band 
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data alone has been tested with satisfactory results ([Chelouche 


& Daniel||2Q12t [Chelouche et al.||2012t l Edri et al.||2012| |Pozo 


Nunez et al.||2013| Chelouche & Zucker||2013| ). In this method, 

two suitable chosen broad band filters are used to trace the con¬ 
tinuum variations and to catch the emission line and continuum 
with the removal of the continuum performed in the cross corre¬ 
lation domain. 

PGC50427 has been classified as a Seyfert 1 galaxy and 
is located at a distance of 102 Mpc ( [Veron-Cetty & Veron| 
|2010| ). In this paper we present the results of a multi-year mon¬ 
itoring campaign carried out on the nucleus of the Seyfert 1 
galaxy PGC50427. We use Photometric Reverberation Mapping 
in combination with dust-reverberation mapping to determine 
the black hole mass, the size of the BLR and dust torus for the 
first time in this source. 


2. Observations and data reduction 

The optical and near-infrared monitoring campaign of 
PGC50427 were carried out using the VYSOS-6, BEST-II, 
BMT and IRIS telescopes (for more details see Sections 2.1-2.2) 
located at the Universitatssternwarte Bochum observatory, near 
Cerro Armazones, the future location of the ESO Extreme Large 
Telescope (ELT) in Chil^ 

In addition to the photometric observations, one single epoch 
spectrum was acquired using the Robert Stobie Spectrograph 
(RSS) at the Southern African Large Telescope (SALT). The 
timeline of all the observations in our campaign is shown in 

Fig[T] 

PGC50427 lies at redshift z = 0.0236, therefore the Ha emis¬ 
sion line falls into the SII 6721 + 30 Aand 6700 + 60 Anarrow- 
band (NB) filters. The characteristics of the source and the galac¬ 
tic foreground extinction values are listed in Table 1. Eigure 2 
shows the position of the NB filters with respect to the Ha emis¬ 
sion line together with the effective transmission of the other op¬ 
tical filters used. 

2.1. Optical monitoring 


VYSOS-6 

Broad-band Johnson B (4330 A), Sloan-band r (6230 A), 
and narrow-band SII (6721 + 30 A , the position of the redshifted 
Ha line) observations were carried out during a monitoring 
campaign between Eebruary 18 and September 01 of 2011, 
with the robotic VYSOS-6 telescope. The VYSOS-6 telescope 
consist of two 15 cm refractors (Takahashi TO A 150E Ortho- 
Apo-chromat Triplet) installed on a common equatorial mount 
(Bisque Paramount ME). Both refractors are equipped with an 
Apogee ALTA U16M 4096 x 4096 pixel CCD, providing a field 
of view (EoV) of about 2.7° x 2.7°. More information about the 
telescope and the instrument has been published by |Haas et al. 

( [20T2I ). 

The data reduction was standardized, including bias, dark 
current, fiatfield, astrometry and astrometric distortion cor¬ 
rections performed with IRAI0 in combination with SCAMP 

^ http://www. astro.ruhr-uni-bochum. de/astro/oca/ 

^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As¬ 
tronomy (AURA) under cooperative agreement with the National Sci¬ 
ence Foundation. 


( |Bertin| |7006| ) and SWARP ( |Bertin et al.||2002} routines. A 
7'.'5 diameter aperture was used to extract the flux, which was 
normalized and converted into photometric fluxes through the 
comparison with 17 non- variable stars located on the same 
images within 30' around the AGN and of similar brightness as 
the AGN. Absolute calibration was performed using standard 
reference stars from |Landolt| ( |2009| ) observed o n the same 
nights as the AGN, considering the at mospheric (|Patat et al 


20111 ) and galactic foreground extinction ( [Schlafly & Finkbeiner 


201 1| corrections. The reader is referred to |Pozo Nunez et ^ 


( |2013| ) for further information about reduction and light curve 
extraction process. 


BEST-II 

Broad-band B (4330 A) and R (7000 A) Johnson observa¬ 
tions were carried out during a monitoring campaign between 
May 01 and July 30 of 2013, and between Eebruary 03 and Au¬ 
gust 19 of 2014, with the robotic 25 cm Berlin Exoplanet Search 
Telescope-II (BEST-II). The BEST-II telescope is equipped with 
a Peltier- cooled 4096 x 4096 pixel Einger Lakes Imager CCD 
KAE-16801, yielding a field-of-view of 1.7° x 1.7° with a pixel 
size of 9 fim. More information about the telescope and the 
instrument has been published by |Kabath et al.| ( [2009] ). The data 
were reduced and calibrated following the same procedures as 
for the VYSOS-6 telescope. 


BMT 


Narrow-band (6700 + 120 A, hereafter denoted as 670) Ashahi 
observations were carried out simultaneously to the broad-band 
monitoring campaign in 2014, using the robotic 40 cm Bochum 
Monitoring Telescope (BMT). The BMT telescope is equipped 
with a two-stage thermoelectric cooled 3072 x 2038 pixel CCD 
SBIG STL-6303, yielding a field-of-view of 41' x 27' with 
a pixel size of 9 pm. The data were reduced and calibrated 
following the same procedures as in the broad-band analysis. 
More information about the telescope and the instrument has 
been published by |Ramolla et al.|p013| ). 


SALT 


The optical spectrum of PGC50427 was observed using 
the 11m Southern African Large Telescope (SALT) on May 10 
2013 (Proposal Code: 2013-1-HETGU-001). The observations 
were performed with the Robert Stobie Spectrograph (RSS) 
mounted at the prime focus of the SALT. We used the PG0900 
grating with a resolving power of 1065 at 6050 A, and the 
spectrum covers a wavelength range between 3200 - 9000 A. 
The spectrum was taken in two identical consecutive exposures 
of 10 min each through the 2" x 8' longslit PL 0200N001 
at a parallactic angle. The detector was operated in normal 
readout mode with a 2 x 2 binning. A Xe spectrum was obtained 
after the object exposure for wavelength calibration. Eor flux 
calibration we observed the standard star G24-9. After the bias 
subtracted file provided by the SALT pipeline, we used standard 
IRAE routines for flat-field correction, cosmic ray rejection, 
2D- wavelength calibration, night sky subtraction and flux 
calibration. Eor the 1-D spectrum, we combined 7 columns, 
corresponding to 1.7738" (0.1267" per unbinned pixel). The 
reduced spectrum is shown in Eig[^ 
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Fig. 2. SALT spectrum of PGC50427. For illustration, the band passes of the filters used for the photometric monitoring are shown (blue 5-band, 
green 5-band, orange r^-band, red 670-band, and brown for Sll-band). The NB 670 and SII catches the redshifted Her line, its fiux is composed by 
the contribution of about 85% Her line and 15% continuum. Note that for actual fiux calculations the filter curves are convolved with the quantum 
efficiency of the Alta U16 and SBIG STL CCDs cameras. 




Fig. 3. Observed light curves of PGC50427, as well as for some of the reference stars in the field of view for the period between February 2011 
and September 2011 (left) and for the period between February 2014 and August 2014 (right). The light curves are vertically shifted for clarity. 
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2.2. Near-infrared monitoring 


IRIS 

Near-infrared (NIR) J (1.25 //m) and Ks (2.15 fim, here¬ 
after denoted as K) observations were carried out between May 
20 and August 11 of 2013 using the 0.8 m Infrared Imaging 
System (IRIS) telescope. IRIS is equipped with a HAWAII-1 
nitrogen-cooled detector array with 1024 x 1024 pixels, yielding 
a field-of-view of 12.5' x 12.5' and a resolution of 0.74'.Vpixel 
( [Hodapp et al.|[MT0| ). Images were obtained by combining 
double cycles of 20 seconds exposure time acquired with 
the observing sequence object-sky-object. The images were 
reduced using IRAF routines. Because the sky background 
emission contribution is one of the most difficult step in NIR 
data reduction, the sky frames observed close to the AGN 
were subtracted from each science frame before ffat-ffeld 
and further corrections. One ffnal image, resulting from the 
combination of all individual frames, is obtained in order to 
remove cosmic rays, hot pixeles and negative residuals from the 
sky-subtracted science frames. The data reduction steps after 
the sky background subtraction and correction for cosmic rays, 
hot pixels and negative residuals, are the same as outlined for 
the VYSOS-6 telescope, including astrometry and astrometric 
distortion correction with SCAMP and SWARP. Light curves 
were calculated relative to 6 non-variable stars located on the 
same ffeld having similar brightness as the AGN. Photometric 
calibration was achieved by using 4 high-quality ffag (AAA) 
Two Micron All Sky Survey (2MASS) stars appearing in 
the same ffeld as the AGN. As already noted for the optical 
treatment of the data ( |Pozo Nunez et al^|2013| ), analysis for 
different aperture photometry was performed considering the 
proper minimization of the host-galaxy contribution and a 770 
diameter aperture was chosen for further analysis. 


3. Results and discussion 

3.1. Optical light curves and BLR geometry 

The optical light curves for campaigns 2011 and 2014 are shown 
in Fig[^ The light curves are published at the CDS. In the cam¬ 
paign of 2011, the 5-band shows a gradual ffux increase by 20% 
from the beginning of March until a maximum is reached at the 
beginning of April (Fig.3, left). After this overall maximum, the 
fluxes undergo an abrupt drop by about 40% until the end of 
May. Between the beginning of June and the end of August 2011, 
the ffuxes are steady. The r-band light curve follow the same fea¬ 
tures than the 5-band light curve albeit with a smaller variability 
amplitude. The latter is expected due to the larger constant host- 
galaxy contribution. Moreover, the r-band also contains a con¬ 
tribution from the strong Ha emission line (Fig[^. The narrow 
Sll-band light curve follows the continuum dominated broad¬ 
band light curves qualitatively, but we can see a 20 day delay 
compared to the continuum variability at the minimum observed 
at the end of June. 

As already discussed in previous PRM studies, the narrow- 
band contains, in addition to the Ha line, a contribution of the 
varying AGN continuum, which must be removed before apply- 
ing cross correlation techniques ([Haas et al.|20TH|Pozo Nunez| 
|et al.|201^ |Pozo Nunez et al.|2013| ). In order to determine this 

contribution, we used the SII and r-band ffuxes, previously cali¬ 
brated to mJy, as is shown with the ffux-ffux diagram in Figure 4. 



rs flux [ mJy ] 

Fig. 4. Flux-flux diagram for the SII and r filter. Black dots denote the 
measurement pair of each night during campaign 2011. The red and 
green lines represent the average flux in the SII and r band respectively. 
Fluxes were measured using circular 7.5" apertures. The data are as 
observed and not corrected for extinction. 



Fig. 5. Same as Fig|^ but for 670 and R data obtained during campaign 
2014. 


The Ha line is contributing, on average, about 70% of the total 
ffux enclose in the Sll-band, while the continuum contribution 
(r-band) is about 30%. Following the usual practice of narrow- 
band PRM, we construct a synthetic Ha light curve by subtract¬ 
ing a fraction of the r-band light curve (Her = SII - 0.3 r), as 
ihustred in Fig[^ The Her light curve was used afterwards to esti¬ 
mate the time delay. For this purpose, we used the discrete corre¬ 
lation function (DCF, [Edelson & Krohk|1988] ) to cross-correlate 
the continuum and the synthetic Her emission line, taking into 

Article number, page 5 ofpff] 


















A&A proofs: manuscript no. dust_BLR 



Fig. 6. Her light curve (red) obtained after the continuum subtraction 
process performed on the SII light curve (black) during campaign 2011. 
The continuum fraction was obtained from the flux-flux diagnostic on 
the SII and r bands; see text for details. 



Fig. 7. Same as Fig|^ but for 670 and R data obtained during campaign 
2014. 


account possible bin size dependenc}r](|Rodriguez-Pascual et al. 

|T^ . 

The centroid of the cross-correlation between the 5-band 
and Ha light curves shows a time delay of 20.5 days (Fig 8). 
Uncertainties in the time delay were calculated using the flux 

^ We note that the choice of a lower or higher time-bin size does not 

change the results if a well sampled data is used dPozo Nunez et al.l 

- 
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Fig. 8. Cross correlation of B and Her light curves for campaign 2011. 
The dotted lines indicate the error range (+lcr) around the cross cor¬ 
relation. The centroid was calculated above the correlation level at 
r > O.^rfnax- The histogram shows the distribution of the centroid lag ob¬ 
tained by cross correlating 2000 flux randomized and randomly selected 
subset light curves (FR/RSS method). The black shaded area marks the 
68% confldence range used to calculate the errors of the centroid. 
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Fig. 9. Same as Fig[^ but for campaign 2014. 


randomization and ra ndom subset selection method (FR/RSS, 
Peterson et al.||2004] ). From the observed light curves we cre¬ 
ate 2000 randomly selected subset light curves, each containing 
63% of the original data points (the other fraction of points are 
unselected according to Poisson probability). The flux value of 
each data point was randomly altered consistent with its (normal- 
distributed) measurement error. We calculated the DCF for the 
2000 pairs of subset light curves and the corresponding centroid. 
From this cross-correlation error analysis, we measure a median 
lag of Tcent = 20.4 ^ BfHa. Correcting for the time dilation fac¬ 

tor (1 -I- z = 1.0236) we obtain a rest frame lag of 19.9 + 0.68 
days BfHa. 

In the campaign of 2014, we see a steep 5-band rise at the 
beginning of April by about 20% which is also observed in the 
5-band but with a smaller amplitude (Fig.3, right). In contrast to 
the steep 5-band flux increase, the narrow 670-band rise appears 
stretched until the end of April with an amplitude of about 13%, 
suggesting a time delay of the Ha line of 15-20 days. After a pe¬ 
riod of constant flux, the 5 and 5 light curves show a small sharp 
bump in June of about 10%, which is observed about 15 days 
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Fig. 10. SALT spectrum of PGC50427, zoomed onto the Her line. The 
dotted black curve represents the spectrum after subtracting the narrow 
[Nil]4/16548,6583 emission lines. The narrow Her and [Nil] models are 
shown in dotted lines at the bottom of the Her profile. The NB 670 and 
SII filters efectively covers the line between velocities -3478km/s and 
+ 1881km/s (red dotted lines) and -1201km/s and +1479km/s (blue dot¬ 
ted lines) respectively. 
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Fig. 11, R BLR - L relationship from data of Bentz et al. ( 2013| ) (black 
dots) with a fitted slope a = 0.533 (dotted line). The diagram also con¬ 
tain the objects from previous photometric reverberation mapping cam¬ 
paigns { [Haas et al.|2011| green dots; |Pozo Nunez et al.|2012f 2013, blue 
stars, and this article, red filled circle). 


later in the 670-band. However, due to insufficient data obtained 
at this period in the 670-band, the delay cannot be measured as 
precise as for the first. We have isolated the Ha emission line by 
the subtraction of the underlying continuum following a similar 
procedure as for the 2011 campaign. The Her line is contributing, 
on average, about 60% of the total flux enclosed in the 670-band, 
while the continuum contribution (/?-band) is about 40% (Fig 5). 
We construct a synthetic Her light curve by subtracting a frac¬ 
tion of the /?-band light curve (Ha = 670 - 0.4 R), as ihustred 
in Fig 7. The centroid from the cross-correlation between the B- 
band and Her shows a time delay of 19.3 days (Fig 9). From the 
FR/RSS cross-correlation error analysis, we measure a median 
lag of Tcent = 18.7 Correcting for the time dilation factor we 
obtain a rest frame lag of 18.3 + 1.03 days. The deduced values 
for the time delay obtained at two different epochs are in quali¬ 
tative agreement, and consistent with the small luminosity varia¬ 
tions in the AGN (~ 10%) accross the monitoring campaign (see 
Section 3.3). We note that the relatively prominent contribution 
of the Hy emission line (Fig. 2) to the 5-band (~ 10%) may re¬ 
sult in a slightly higher time lag centroid from the DCF as the B- 
band light curve contains not only the AGN continuum but also 
the lagged Hy signal. In the net effect, our estimate B/Ha may 
underestimate the true B/Ha time delay. However, the current 
data does not allow to quantify this effect; this would requiere 
the use of emission line free continuum bands, for instance with 
the ultraviolet band. 


3.2. Central black hole mass 

Assuming that the BLR emiting gas clouds are in virialized mo¬ 
tion around the central black hole, the mass of the black hole 
can be estimated as Mbh = / • ^blr • o-y/G. The velocity cry of 
the emission-line region is determined from the line dispersion 
(t^iine) or from the full width at half maximum (FWHM) of the 
line profile. The scaling factor / depends on the geometry and 
kinematics of the BLR. Most of the results presented in previ¬ 
ous reverberation studies have been carried out considering only 
the virial product cro-y/G, ie. assuming a scaling factor f = I 
(Peterson et al. 2004, and references therein). 


The broad Ha emission line is blended with the narrow Ha 
and [NII]/1/16548,6583 narrow emission lines. In order to re¬ 
move the narrow components, we model the [SII]/1/16716,6731 
double t with a multi-Gaussian profile as described in 
|& Ho| ( |2004| ). The model is shifted and scaled to fit 
-f [Nil]446548,6583 narrow lines and subtracted from the ob¬ 
served broad Ha line profile as shown in Fig. 10. The ratio of 
the [Nil] lines is fixed at the theoretical value of 2.96 and the 
relative positions of the narrow Ha and [Nil] lines are deter¬ 
mined by their laboratory wavelenghts. After removing the nar¬ 
row emission lines, the Ha profile was isolated by the subtraction 
of a linear continuum fit, obtained through interpolation between 
two continuum segments on either end of the line. Figure 10 il¬ 
lustrate the original Ha emission line profile together with the 
subtracted narrow emission line profiles. The velocity dispersion 
after removal of narrow lines is crune = 1020 km s“\ which has 
been corrected for instrumental velocity dispersion to obtain an 
intrisic profile width. The feasibility of the use of single-epoch 
(SE) spectra for the black hole mass determination has been es¬ 
tablished in previous investigations. On average, uncertainties 
of -30% have been reported for black hole mass determination 
from single epoch spectra measurem ents (e.g. |Vestergaard|20'02| 
Woo et ffi.|2007||Denney et al.|2009| ). 

Using the derived time delay r = 18.3d for epoch 2014 and 
the velocity dispersion (with 30% uncertainty), the virial black 
hole mass is M^triai = (3 ±2) x lO^M©. Considering the factor 
/ = 5.5 + 1.8, based on the asumption that AGNs follow the 
same Mbh - cr* relationship (Onken et al. 2004), we determine 
a central black hole mass Mbh = (17 + 11) x 10^ M©. 

Assuming a symmetric BLR, the dimensionless factor / de¬ 
pends on the unknown inclination of the BLR (/ = 

The optical and the Ha emission line variability observed in 
PGC50427 suggest a disk- like BLR geometry with low incli¬ 
nation i < 30° (Pozo Nunez et al., in preparation), hence the 
geometry-scaling factor / may be much higher than the com¬ 
monly used f = 5.5. Therefore the black hole mass derived here 
should be considered as lower limit. 


Greene 
the Ha 
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Fig. 12. Flux variation gradient diagrams for 2011 (left), 2013 {middle) and 2014 {right) campaigns. The solid li nes delineate the bise ctor regression 
model yielding the range of the AGN slope. The dashed lines indicate the range of host slopes determined by [Sakata et al.| ( [20T0] ) for 11 nearby 
AGN. The intersection between the host galaxy and AGN slope (red area) gives the host galaxy flux at the time of the campaign in both bands. 


3.3. Host-subtracted nuclear luminosity and the BLR size - 
luminosity relationship 

To determine the AGN luminosity free of host galaxy contribu¬ 
tions, we applied th e flux variation gradie nt (FVG) method, orig¬ 
inally propose d by |Choloniewski| ( 1981| ) and later modified by 
[Winkler et al.| ( 1992) . A detai led description of the FVG method 
on PRM data is presented in [Pozo Nunez et'ST (2012). In this 
method the fiuxes obtained through different filters and same 
apertures are plotted in a fiux-fiux diagram. The fiuxes follow 
a linear slope representing the AGN color, while the slope of 
the nuclear host galaxy contribution (including the contribution 
from the narrow line region (NLR)) lies in a well defined range 
(0.4 < < 0.53, for 8'.'3 aperture and redshift z < 0.03, 

Sakata et al.[[2010| ). The AGN slope is determined through a 


linear regression analysis. Averaging over the intersection area 
between the AGN and the host galaxy slopes yields the actual 
host galaxy contribution at the time of the monitoring campaign. 
Figure shows the FVG diagram for the B and r fiuxes cor¬ 
responding to campaign 2011 and for the B and R fiuxes cor¬ 
responding to campaign 2013/2014 obtained during the same 
nights and through a 7'.'5 aperture. All the fiuxes has been cor¬ 
rected for galactic foreground extinction. 

The bisector linear regression method yields a linear gradient 
of Tbr, = 1.18 + 0.06, Tbr = 1.22 + 0.03 andF^^ = 1.10 + 0.05 
for campaigns in 2011, 2013 and 2014 respectively. The results 
are consistent, within the uncertainties, with the gradients ob- 
tained for other Sey fert 1 galaxies by [Winkler et al.[ ( [1992| ) and 
[Sakata et ar| ( [2010[ ). The AGN fiuxes at the time of the monitor¬ 
ing can be determined by subtracting the host galaxy fiuxes from 
the total fiuxes. The host galaxy subtracted average AGN fiuxes 
and the host galaxy flux contribution of PGC50427 are listed 
in Table Also listed in Table are the interpolated rest frame 
5 lOOA fiuxes and the monochromatic AGN luminosity /IL^(agn) 
at 5100A obtained at the distance of 102 Mpc. The rest frame 
flux at 5100A was interpolated from the host-subtracted AGN 
fiuxes in both bands, assuming for the interpolation that the AGN 
spectral energy distribution (SED) is a power law {Fy oc y^) with 
a spectral index a = log{fxAGN/fyAGN)/ log{yx/yy), where 
and yy are the effective frequencies in the X and Y bands, re¬ 
spectively. The error was determined by interpolation between 
the ranges of the AGN fiuxes ±cr in both filters, respectively. 

The position of PGC50427 on the BLR size-luminosity di¬ 
agram is shown in Figur e pG] The values of other galaxies are 
taken from [Bentz et aL] ( [2013[ ) and from previous photometric 
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Fig. 13. Optical and NIR host galaxy subtracted light curves obtained 
between May and August 2013. 


reverberati on mapping campaigns ( [Haas et al.[201 l[[Pozo Nunez 
et al.|2012] 2013). For this figure we converted the measured Htr 
size into the size of the HjS BLR using the weighted mean ra¬ 
tio for the time lag T{Ha) : t{HI5) : 1.54 : 1.00, obtained by 
Bentz et al. (2010) from the Lick AGN Monitoring Program of 
11 low-luminosity AGN. 


3.4. Infrared light curves and Dust-torus size 

Figure [^depicts the optical and near-infrared normalized light 
curves of the nucleus of PGC50427 obtained during the 2013 
campaign. The light curves are published at the CDS. To decon¬ 
volve the host galaxy and the nuclear flux contributions, we used 
the flux v ariation gradient (FVG) m ethod in the same way as de¬ 
scribed in [Pozo Nunez et al.] ( [2014[ ). The spectral energy distribu¬ 
tion of the variable component remain constant with time, and in 
consequence the slopes obtained from the optical and NIR flux 
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Table 2. Total, host galaxy and AGN optical fluxes given in mJy for the different campaigns. 


Campaign 

Filter 

Total 

(mJy) 

Host 

(mJy) 

AGN* 

(mJy) 

fAGN((l + Z)5100A 
(mJy) 

lOOA 
erg s~^) 

2011 

B 

2.71+0.12 

0.91+0.17 

1.80+0.20 

1.66+0.20 

1.18+0.14 


rs 

3.95+0.11 

2.42+0.15 

1.53+0.19 



2013 

B 

2.36+0.10 

0.83+0.15 

1.53+0.18 

1.40+0.17 

1.00+0.12 


R 

3.56+0.11 

2.35+0.12 

1.21+0.16 



2014 

B 

2.48+0.10 

0.80+0.13 

1.68+0.16 

1.72+0.17 

1.22+0.12 


R 

3.94+0.12 

2.12+0.15 

1.82+0.19 




a) 

AGN fluxes values /agn = ftotal - fhost with uncertainty range ctagn = i(^]otai + 




Fig. 15. Distribution of AGN NIR J-K colors of Seyfert galaxies taken 
from |Glass|2004| after host galaxy subtraction. The red vertical solid line 
correspond to the / - ^ color of PGC50427. The filled black squares 
indicate Blackbody color temperatures. 


Fig. 14. Flux variation gradient diagram in the NIR. The data are rep¬ 
resented by the black dots. The solid blue lines represent the bisector 
fit, yielding the range of the AGN slop es. The dashed red lines indicate 
the range of host slopes determined by |Suganuma et al.| ( |2006| ). The in¬ 
tersection between the host galaxy and AGN slope (red area) gives the 
host galaxy flux in the respective bands. 


ratios (B/R and J/K) allow us to separate the AGN flux through 
the use of a well-defined range of host galaxy colors ( |Glass|2004| 
[Sakata et al.||2010| ). Figure shows the J flux plotted against 
the K flux for each night. The bisector linear regresion method 
yields a linear gradient of F//^ = 0.38 + 0.05 which correspond 
to a nuclear color J-K = 1.97 + 0.07, and is consistent with 
the FVG colo rs of Seyfert galaxies {J - K ^ 2.3) measured by 
Glassj ( |2004| ). Figure shows the distribution of / - ^ col¬ 
ors of Seyfert galaxies previously studied by [Glassj p004[ ) and 
the value for PGC50427 obtained in this work. The color of 
PGC50427 correspond to a blackbody temperature of 1850K, 
which is in agreement with the sublimation temperature of hot 
pure graphite grains ( [Mor & Netzer|2012| ), and supports the pre¬ 
vious conclusions that IR emission observed from Seyfert galax¬ 
ies is dominated by thermal radiation from the hot dust close to 
the central AD ( [Suganuma et al.|2006 ). 


The optical B and R-bands, which are mostly dominated by 
the AGN continuunj^ show a similar varibility behavior. The 
NIR J and ^-bands, which are dominated by the innermost hot 
dust reacts with delayed response to the AGN continuum. 

The time delay between the AGN continuum and dust emis¬ 
sion can be estimated by cross-correlation of the optical and NIR 
light curves yielding the average radius of the innermost dust 
torus. We correlated both the B- and ^-band light curves and the 
B- and J- band light curves using the discrete correlation func¬ 
tion (DCF, [Edelson & Krohk[1988| ). The cross correlation of BjJ 
shows two peaks, one small correlation peak around lag 19 days 
obtained above the correlation level at r > OAr^^ax and a major 
peak with a lag of 46 days obtained above the correlation level at 
r > as shown in FigureSimilar features can be seen 

in the cross correlation of BjK, a small correlation peak with a 
lag of 19 days, a major peak with a lag of 47 days, and an ad¬ 
ditional third and smallest peak with a lag of 67 days as shown 
in Figure The NIR host galaxy corrected light curves show 
features of similar amplitude and sharpness as the optical host 
galaxy corrected light curves. In addition, the cross correlation 
are at nearly zero level at lag 0 days. One expects that if the dust 
distribution is spread at different line-of-sight distances, the ob¬ 
served echo will be smeared out in time and the J and K light 
curves will show a smoother variability. Therefore the observed 


^ We note that the R-band also contains a small contribution from the 
strong Her emission line. 
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Fig. 16. Cross correlation of B and J light curves. The dotted lines in¬ 
dicate the error range (+lcr) around the cross correlation. The centroid 
was calculated above the correlation level at r > 0.6rmax- The histogram 
shows the distribution of the centroid lag obtained by cross correlat¬ 
ing 2000 flux randomized and randomly selected subset light curves 
(FR/RSS method). The black shaded area marks the 68% conhdence 
range used to calculate the errors of the centroid. 
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Fig. 18. Lag - luminosity relationship based on the data of 
|et al.| ( [2006] ) and|Koshida et al.|(|2014| ). The solid line is the 
gression from |Koshida et al.| ( |2014| ). PGC50427 (green dot) lies close to 
the regression line. The dashed blue line indicates the dust s ublimation 
radius expected at a given nuclear luminosity My (from|Kishimoto| 
let al.|2007| ). 
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Fig. 17. Same as Fig[^ but for B and K light curves. 


NIR light curves of PGC50427 argues in favor of a face-on torus 
geometry. 

The uncertainty of the lag time r was estimated using the 
flux randomization and random subset-selection method (FR and 
RSS, |Peterson et al.|2004| ). The median of this procedure yields 
Tcent = 46.8 days and Teem = 47.8 ^3 2 days for B/J and B/K, 
respectively. Correcting for the time dilation factor, we obtain a 
rest frame time-delay Trest = 45.7 + 1.47 days and Trest = 46.7 + 
2.15 days for BjJ and BjK, respectively. 

Dust-reverberation studies of Seyfert 1 galaxies have shown 
that the dust torus size obtained by the cross-correlation of the 
optical V and ^-bands is proportional to the square root of the 
optical luminosity r cx (Suganuma et al. 2006, Koshida et al. 
2014). Figure shows the position of PGC50427 on the Dust 
size-luminosity diagram. To obtain the optical V-band fluxes we 
interpolated the AGN B and R band fluxes. 
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4. Summary and conclusions 

We have performed three monitoring campaigns between 2011 
and 2014, as well as obtaining a SALT spectrum, of the Seyfert 1 
galaxy PGC50427. We determined the broad line region size, the 
dust torus size, the virial black hole mass and the host-subtracted 
AGN optical and NIR luminosity. The results are: 

1. The cross-correlation of the Her emission line with the op¬ 
tical continuum during campaign 2011 yields a rest-frame 
time delay Trest = 19.9 +0.68 days. During campaign 2014 
the cross-correlation of the Her emission line with the optical 
continuum yields a rest-frame time delay Trest = 18.3+1.03. 

2. With the velocity dispersion obtained from a single epoch 
spectrum in 2013, we determine the black hole mass Mbh = 
(17 + 11) X 10^ Mq assuming a geometrical factor of 5.5. 

3. Using the flux variation gradient method (FVG) we deter¬ 
mine the host galaxy subtracted optical AGN luminosity of 
PGC50427 at different epochs of observations Lagn- 20 ii = 
(1.18 + 0.14) X \Q^^ergs-\ Lagn -2013 = (1-00 + 0.12) x 
W^ergs-\ and Lagn = ( 1-22 + 0 . 12 ) x \0^^ergs-\ 

4. From the NIR light curves in 2013, we determine a lag time 
of Trest = 45.7 + 1.47 days and Trest = 46.7 + 2.15 days for 
BIJ and BjK, respectively. The relatively sharp dust echo 
observed in the NIR light curves argues in favor of a face- 
on torus geometry. The infrared lag time is correlated with 
the optical luminosity according to Rdust ^ in agreement 
with previous investigations. The infered inner size for the 
dust torus in PGC50427 suggest that the location of the ther¬ 
mal emitting region is located well outside the BLR which 
support the unified scheme of AGNs. 
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